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ABSTRACT 

This paper presents a comparison of the measured 
and the computed spatial stability properites of an 
aspect ratio 2 supersonic shock free elliptic jet. The 
shock free nature of the elliptic jet provides an ideal 
test of the validity of modeling the large scale co- 
herent structures in the initial mixing region of non- 
drcular supersonic jets with linear hydrodynamic sta- 
bility theory. Both aerodynamic and acoustic data 
were measured. The data are used to compute the 
mean velocity profiles and to provide a description of 
the spatial composition of pressure waves in the ellip- 
tic jet. A hybrid numerical scheme is applied to solve 
the Rayleigh problem governing the inviscid linear 
spatial stability of the jet. The measured mean veloc- 
ity profiles are used to provide a qualitative model for 
the cross sectional geometry and the smooth velocity 
profiles used in the stability analysis. Computational 
results are presented for several modes of instability 
at two jet cross sections. The acoustic measurements 
show that a varicose instability is the jet’s perferred 
mode of motion. The stability analysis predicts that 
the Strouhal number varies linearly as a function of 
axial distance in the jet’s initial mixing region, which 
is in good qualitative agreement with pre..ous mea- 
surements. 

1.0 INTRODUCTION 

For years the origin of noise production by super- 
sonic round jets has been pursued vigorously both 
theoretically and experimentally. Several theoretical 
models exist to predict this noise, but the most suc- 
cessful concept to date is the quasi-linear instability 
wave model. Tam, ref. 1, has applied this model 
to provide a comprehensive prediction for both jet 
and shock noise components. A natural extension 
of this model is its application to non-circular super- 


sonic jets. Morris, ref. 2, applied the instability wave 
model to elliptic jets and was able to determine that 
as the exit plane momentum thickness increases, the 
dominant spatial mode of instability occurs along the 
major axis. While this work provided a foundation 
for analysis of non-circular jets, only those velocity 
profiles that allowed a separable solution were appli- 
cable. Thus the analysis could not be performed for 
those flows whose jet column deformed as it evolved 
downstream. To overcome this important limitation, 
Baty and Morris, ref. 3, introduced a scheme to com- 
pute non-separable solutions of the stability equation 
which occur when the mean velocity profile is de- 
formed as the flow evolves downstream. This method 
is based on a conformal mapping of the irregularly 
shaped jet cross section along with a hybrid numeri- 
cal scheme. 

In the present paper the hybrid numerical scheme 
is applied to solve the Rayleigh equation which gov- 
erns the inviscid linear spatial stability of an elliptic 
jet with an initial aspect ratio of 2. An elliptic noz- 
sle was designed for this study to produce a shock 
free flow of Mach 1.5. The spatial modes of instabil- 
ity and neutral instability frequencies as a function of 
axial location were calculated. An experimental pro- 
gram was also conducted to determine the preferred 
spatial mode of instability and identify the Strouhal 
number associated with peak amplitude noise radia- 
tion from the jet. Previous studies have shown that 
the frequency for the most unstable shear layer mode 
is directly related to this Strouhal number. Mean ax- 
ial velocity profiles and contours were also obtained 
to improve the numerical model. Only data associ- 
ated with the operation of the supersonic nozzle at 
its shock free design point are presented in this pa- 
per. 
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The elliptic nozzle tested was designed to produce 
Mach 1.5 axial flow at the exit. The internal geome- 
try is presented in figure 1. The convergence section 
consisted of a circular to elliptic transition leading to 
an elliptic subsonic contraction of aspect ratio 2. This 
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aspect ratio was held constant throughout the nossle 
divergence section to the exit where the dimensions 
are 2.000 by 1.000 inches. The axial length of the 
divergence section is 5.015 inches. 

A nossle design method developed by Seiner, et al., 
ref. 4, using both one and two dimensional inviscid 
isentropic theory, was applied to generate the interior 
wall coordinates. The design goal was to produce uni- 
form Mach 1.5 parallel flow at the nossle exit where 
the pressure waves in the plume differed from the am- 
bient pressure by less than 10 percent. Shock waves 
that satisfy this pressure minimisation criterion pro- 
duce insignificant levels of shock noise relative to other 
jet noise sources. Since the nossle design procedure is 
inviscid, the r ossle was operated at a nossle pressure 
ratio of NPE. = 3.78 to account for boundary layer 
effects. At this NPR the flow only contained weak 
waves with an exit Mach number of 1.52. 

Aerodynamic measurements were performed in the 
NASA Langley Jet Noise Laboratory. These mea- 
surements utilised total and static pressure probes de- 
signed for supersonic flow. The total pressure probe 
consisted of a 0.080 inch diameter tube with a 0.020 
inch sensor port where the probe inlet was machined 
flat to create a well defined detached shock when sub- 
merged in supersonic flow. The static pressure probe, 
based on the design of Pinckney ref. 5, was a dual 
conical type which recovers the tip static pressure 
faster than conventional static probes. This allows 
the holes to be positioned near the probe tip which 
is important for flows with large axial pressure gra- 
dients. Both probes were installed in a probe wing 
designed for supersonic flow which was attached to a 
remote controlled three dimensional traversing mech- 
anism capable of 0.001 inch movements. Centerline 
data were collected as well as one dimensional pres- 
sure proves and two dimensional pressure maps at 
various axin' positions in the jet. 

A schlieren system was assembled in the Jet Noise 
Lab to give an indication of the density gradient field 
in the flow. A vertical knife edge was used to give the 
effect of the density gradient component in the direc- 
tion of the flow. Figures 2a and 2b show the schlieren 
optical records in planes parallel to the major and mi- 
nor axes respectively, ref. 4. Waves associated with 
the weak variatioaM K axial static pressure can be 
seen in the figu i ■ to the high sensitivity of the 
schlieren apparat us. 

Also collected in this facility were near-field acous- 
tic measurements using an elliptic microphone array, 
figure 3. The array consisted on twelve Bruel and 
Kjaer 1/4 inch pressure micorphones positioned to 
form an aspect ratio 2 ellipse whose major axis length 
was 5.0 inches. The positions of the microphones in 


Cartesian coordinates (origin at the jet axis) are given 
in parenthesis in figure 3 for a single quadrant. The 
microphones were positioned perpendicular to the jet 
axis in the nostle exit plane and spaced at equal ellip- 
tic arc lengths of 1.009 inches on the array. The data 
were gathered using a wide band analog magnetic 
tape recorder band limited to 20 kHs. Spectra were 
acquired in real time via a Scientific Atlanta Model 
SD380 Fast Fourier Transform (FFT) analyser. By 
heating the flow electrically, all measurements were 
obtained at a stagnation temperature of 40 degrees 
Celsius. 

Additional acoustic data were collected in the 
NASA Langley Anechoic Noise Facility, ref. 6. 
Twenty six Bruel and Kjaer 1/4 inch microphones 
were positioned on a linear array parallel to the jet 
axis and perpendicular to the nossle’s minor axis, 
i.e., facing the nossle exit’s wide dimension. Figure 
4 shows the seven radial positions used for the array, 
where D t1 is the elliptic nossle’s equivalent diameter 
relative to a round nossle defined by: 

D. f = 2 BVAR (2.1) 

where B is the nossle minor radius and AR is the nos- 
sle aspect ratio. In the present case D tf is equal to 
1.414 inches. The microphones were spaced at equal 
intervals of 3 equivalent diameters. The fourteenth 
microphone from the top of the array was located 
in the nossle exit plane. The acoustic signals were 
tape recorded using intermedidate band recorders lim- 
ited to 20 kHs and were post-test processed using the 
SD380 FFT analyser. Because no electric heat was 
available in this facility, the stagnation temperature 
varied from approximately 26 to 31 degrees Celsius 
throughout the experiment. 

3.0 EXPERIMENTAL RESULTS 

This section presents both near field acoustic and 
plume mean flow data. Even though measurements 
were obtained for a wide range of nossle pressure ra- 
tios, only the fully expanded results are presented 
here. The linear inviscid Rayleigh equation, as used 
in the present analysis, is not valid for plumes that 
contain shocks. Thus all data presented was obtained 
by operating the elliptic nossle at a nozzle pressure 
ratio of 3.76, which corresponds to a fully expanded 
Mach number of 1.52. 

3.1 Acoustic Measurements 

As described in section 2.0, near field acoustic data 
was obtained using a 26 element linear microphone r»r- 
ray and a 12 microphone element array with sensors 
positioned along an elliptic curve. The linear array 
provided data that enabled determination of impor- 
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tent spectral and directional characteristics of the ra- 
diated sound field. The 12 element array, which had 
sensors distributed equidistant along an ellipse sur- 
rounding the nossle enabled determination of the jet’s 
preferred spatial structure. 

The linear array was positioned at seven radial lo- 
cations from the minor axis, thus providing a total of 
182 measured points. The radial locations varied from 
10 to 40 in S D, t increments. Figure 5 shows con- 
tours of equal overall sound pressure level (OASPL) 
obtained by applying the DI-3000 contour routine to 
the measured data test matrix. The scaled position of 
the nossle is also shown in the figure. The contours 
are presented in 2 dB increments and range from a 
value of 104 dB to 132 dB. 

The contour data of figure S clearly shows that most 
noise emitted from the elliptic nossle’s plume is di- 
rected at a narrow angle to the direction of the flow. 
The majority of noise is contained within a 40 degree 
cone centered at the nossle exit. This noise is of prime 
importance to the analysis contained in this paper. 
The data suggests that this noise originates from a 
region between 5 and 6 diameters from the nossle exit 
and is directed at 30 degrees to the jet axis. The con- 
tour data shows no evidence of screech or broadband 
shock associated noise, thus indicating that the nossle 
design procedure, ref. 4, was successful. The contour 
data results are also consistent with those measured 
by Yu and Dosanjh, ref. 7, for the fully expanded 
Mach 1.5 axisymmetric nossle and those predicted by 
Tam and Burton, ref. 8. One significant difference 
between the round and elliptic jets is that acoustic 
radiation from the elliptic jet has an apparent origin 
closer to the nossle exit than is found for the round 
jet. Based on velocity measurements shown in the 
next section, the convection velocity of the large scale 
turbulent structures in the jet shear layer is subsonic. 
Thus the sound radiation would not include noise due 
to the eddy Mach wave mechanism. 

A typical narrow band spectrum of the noise ra- 
diated into the cone is given in figure 6. The loca- 
tion of the microphone for this spectrum is shown 
in figure 5 as the solid circle in the upper left-hand 
quadrant of the contour. The frequency correspond- 
ing to the spectrum peak amplitude is labelled as /, 
in figure 6. Since the spectrum peak is broad, the 
value f t = 2.4 kHs only represents a nominal value. 
Based on this frequency, however, the nondimensional 
Stronhal number, 

St t = = 0.18 (3.1) 

compares well with tb-. value 51, = 0.22 previously 
reported by Seinr ref. 9, for the Mach 1.5 

axisymmetric jet. 


3.2 Measnred Modes of Instability 

An example nossle exit microphone spectrum is 
shown in figure 7. This spectrum is from the micro- 
phone located along the minor axis in the direction of 
the linear array microphones. Even though the spec- 
trum appears irregular with multiple broadband spec- 
trum peaks, the Strouhal number St c — 0.18 compo- 
nent, identified in figure 6, is clearly present. This 
suggests that the major noise component may also 
emit energy directed upstream from the region be- 
tween 5 and 6 D, n . 

The nossle exit array data was used to evaluate 
the preferred spatial structure of the large scale in- 
stabilities in the jet plume. This was accomplished by 
Fourier decomposion of the asimuthal variation of the 
measured cross correlation coefficients between the 12 
array microphones for several octave bands of noise. 
Figures 8a, 8b, and 8c show the results of this anal- 
ysis for the Strouhal number bands Stc = 0.08, 0.13, 
and 0.16. The cross correlation coefficient amplitudes 
are presented in terms of the normalised arc distance 
around the array. The arc distance is the fundamen- 
tal period for the Fourier decomposition, each micro- 
phone being spaced in equal increments of sue length. 
Microphone 4, which lies along the minor axis is used 
as the reference microphone in computing the cross 
correlation coefficient. Each figure contains summary 
findings of the modal content associated with each oc- 
tave band of noise. Here Ao denotes the varicose or 
axisymmetric mode. The amplitudes associated with 
even symmetry are given by An, and those with odd 
symmetry by Bn. Thus an equal combination of A1 
and Bl modal amplitudes will produce a flapping mo- 
tion in the jet column. 

The modal decomposition presented in figures 8 in- 
dicates that only the varicose and flapping modes have 
sufficient amplitudes to be an interest. All higher or- 
der modes are insignificant. The varicose instability 
surfaces as the dominant mode for all three octave 
bands. The predominance of the varicose mode in- 
creases with increasing Strouhal number, as shown 
by comparing figures 8a, 8b, and 8c. Thus for the 
principle noise producing mechanism, the dominant 
mode of instability is the varicose mode. 

An important point should be made concerning 
Fourier decomposition of the array element cross cor- 
relation coefficients. The Fourier modal amplitudes 
depend on the coordinate surface on which the mi- 
crophones are located. In the present analysis, micro- 
phones are located around the elliptic nossle along an 
elliptic cylindrical coordinate, (ECC). For round jets 
microphones are generally distributed along a circu- 
lar path. The choice of microphone orientation is dic- 
tated by a predetermined understanding of the nature 
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of the noise generation process. For round jets we ex- 
pect the noise source mechanism to be axisymmetric. 
For elliptic jets this type of thinking can be primitive. 
If the jet column remains undistorted, that is, it has 
a spread rate that does not depend on asimuth, then 
the ECC may represent a prudent choice. If however, 
significant distortion of the jet column takes place, 
then another coordinate system needs to found for 
synthesising the fundamental modes. Since only two 
modes were found to contain significant energy, the 
present analysis implies that the choice of the present 
coordinate system may have been correct. In the fol- 
lowing section the flow measurements will also show, 
at least to the end of the potential core, the jet column 
undergoes only minor distortion. 

3.3 Flow Field Measurements 

As discussed above, application of the linear invis- 
dd instability wave model requires empirical input of 
the mean velocity distribution throughout the jet. To- 
ward this end, the mean axial velocity was measured. 
Both *he minor and major axis surveys were made at 
seven udal locations ranging from 0.2 to 10.0 D, t . In 
addition, a high resolution traverse of the centerline 
distribution of axial velocity was made from the nos- 
sle exit to 28 £> M along with four 351 point quadrant 
flow maps at the 0.25, 0.50, 0.75, and 1.00 potential 
core lengths. 

The centerline velodty distribution ... measured to 
determine the length of the jet’s potential core and for 
supersonic jets to gage the adequacy of the nossle de- 
sign method to achieve shock free flow. Figure 9 shows 
the measured centerline velodty when operating the 
elliptic nossle at its design point. The axial variations 
in mean flow velocity in the first 7 D tt are due to im- 
perfections in the nossle design. These variations are 
often stronger with round jets, due to their focusing 
of waves along the jet centerline. Based on this data 
it is possible to determine the axial extent of the po- 
tential core, L c , as that point where the mean velodty 
is 99 percent of the jet exit velodty. This point oc- 
curs at X c — 5.06 D; as is more clearly shown in the 
expanded view of figure 10. 

Figures 11a and 11b show, respectively, 351 point 
axial mean flow maps for X — 0.25 L t and 0.50L t . 
The contours range from a value of 100 to 1400 ft/sec 
in 100 ft/sec increments. In genera] the flow appears 
to distort gradually. 

The initial distribution of momentum thickness has 
been used in the past, for example by Ho and Gut- 
mark, ref. 10, and by Bridges and Hussain, ref. 11, 
as a length scale for the initial turbulent shear layer. 
The asimuthal distribution of the momentum thick- 
ness around the periphery of the nossle exit controls 


the initial role-up of shed vorticity, and the eventual 
distortion of the jet column. The axial distribution of 
the momentum thickness along both the major and 
minor axis then would provide an indication of dis- 
tortion of the jet column with downstream distance. 
For compressible flows the momentum thickness is de- 
fined by 

r* pV V 

'=L w)} l -v. )ir M 

where V 9 is the value of centerline velo ity at a given 
axial location. 

The mean velocity profiles along the major and 
minor axis were integrated numerically for momen- 
tum thickness at the seven axial stations. Figure 
12 shows the resulting axial momentum thickness for 
both major and minor axes. As can be observed, the 
momentum thickness is nearly equal on both major 
and minor axes over the length of the potential cere. 
This suggests that the initial strength of shed vor- 
iicity is constant around the periphery of the nossle, 
thus producing little tendency for distortion of the 
jet column as the flow evolves. In the potential core 
region the momentum thickness also grows linearly, 
a result consistent with the well known behavior of 
round jets. Near the end of the potential core the 
momentum thickness along the major axis grows at 
a much greater rate than that along the minor axis. 
Based on figure 12, this occurs just before the end of 
the potential core. 

4.0 NUMERICAL MODEL 

Consider the shock free jet issuing from the aspect 
ratio 2 elliptic nozsle described in section 2.0, with 
a uniform exit Mach number of 1.52. To establish 
the governing stability equation, assume a Cartesian 
coordinate system with the jet axis aligned in the z 
direction. Let W(z,y) denote the mean axial velocity 
of the jet. The local parallel flow assumption of hy- 
drodynamic stability is applied neglecting the mean 
velocity components in the z and the y directions. 
If this assumption is not made and the effects of flow 
divergence are considered, a problem results which re- 
quires a multiple-scale analysis. Such effects are not 
addressed in the present study. 

A linear, elliptic, partial differential equation for 
the pressure associated with the large scale coherent 
structures is obtained by taking the divergence of the 
momentum equation and using the ideal gas law with 
the equations of continuity and energy. The resulting 
equation is linearised about the mean flow. The veloc- 
ity fluctuations are eliminated in favor of the pressure 
fluctuation using the linearised momentum equations. 
Then assuming that the pressure fluctuation may be 
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written in the form: 

p(*,v,*,i) = p(*,v)** P»7 (41) 

where 

7 = at - ut (4.2) 

a it the axial wave number, and u it the instability 
wave frequency, the equation for p may be expressed 


(A - a .’)p + VK- = 0 (4.3) 

Equation (4.3) is the Rayleigh equation governing the 
in viscid, compressible, isothermal spatial stability of 
the elliptic jet. In (4.3), the compressible wavenum- 
ber, a e , is defined as: 

a e 3 = o* - Mj*(v - aW) 3 (4.4) 

where Mj is the exit Mach number of the jet, which 
is 1.52 in the present study. The equations exhib- 
ited have been non-dimensionalised with respect to 
the pressure, temperature, acd velocity at the jet exit. 
The length scale is based on an equivalent radius de- 
fined by: 

L = y/AB = D.,/2 (4.5) 

where A and B are the major and minor dimensions of 
the jet, respectively. In order to determine the pres- 
sure p, in equation (4.3), boundary conditions must 
be added. In this case p satisfies: 

p — » 0 as |x| — * +oo (4-6) 

and 

p is finite as |x| — » 0 (4.7) 

The computational goal in solving the stability prob- 
lem defined by equations (4.3) to (4.6) is to determine 
the complex wavenumber spectrum, a c , for a given 
frequency, u. 

4.1 Jet Model 

To analyse the elliptic jet it is advantageous to in- 
troduce a new coordinate system. Transformations 
are used to map standard computational rectangles 
onto regions modeling the elliptic jet’s shear layer 
cross sections. Conformal maps are applied to gen- 
erate the necessary transformations. These maps are 
desirable analytically, since they generate diagonal 
metric terms which simplify the computational form 
of the Rayleigh equation. 

Conformal maps may be constructed very efficiently 
for a wide range of simply and doubly connected re- 
gions bounded by elliptic curves, Baty, ref. 12. In the 


present study, the conformal maps will be defined by 
the complex cosine: 

w = acoss for s = y 1 + »y* (4.8) 

where a is a real number. Here equation (4.8) is equiv- 

alent to introducing an elliptic cylindrical coordinate 
system. L* this coordinate system, the Rayleigh equa- 
tion becomes: 

= 0 «- 9 > 

where metric term g is defined by : 

p = |asinx| J (4-10) 

The ECC generate confocal elli~-ic regions which 
capture the qualitative shape of the elliptic jet’s cross 
section. The simple conformal maps defined by (4.8) 
were applied to model the jet’s shear layer geometry 
rather than the more realistic models developed in ref. 
12 for two reasons. Firstly, the use of confocal coor- 
dinate allowed the results of the numerical scheme to 
be validated for an incompressible case given in Mor- 
ris, ref. 2. Secondly, a coordinate system was needed 
which modeled the gross physical characteristics of 
the elliptic jet’s cross sectional geometry. Using the 
ECC the half velocity curve in the shear layer is ap- 
proximated closely. Figure 13 shows a typical jet cross 
section and the half velocity curve in nondimensional 
coordinates. 

Furthermore, it is not clear that using a different 
conformal map would yield a more accurate shear 
layer model. The measured jet velocity contours of 
figures 11 indicate that the'noxsle has a slight imper- 
fection along the wall joint on the minor axis near the 
exit plane. This imperfection was caused by a mis- 
match of the nossle halves. The measured velocity 
profiles show that the nossle distorts the elliptic jet by 
forcing the jet to bulge out on the top half of the jet. 
This phenomenon is shown in figure 11a. The exist- 
ing conformal mapping techniques, ref. 12, developed 
to mode] arbitrary elliptic jets, are unable to gener- 
ate the coordinate transformations needed to analyse 
such small scale flow imperfections. Therefore, the 
most simple coordinate transformations which cap- 
ture the global features of the jet’s cross sections were 
introduced into the numerical model. 

The mean velocity profile used in the computation 
is based on a generalisation of the profile given by 
Michalke, ref. 13. Details of the development of the 
velocity profile in ECC can be found in Morris, ref. 
14. For all the computations the mean axial velocity 
profile is assumed to be given by: 
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W(y>)=i(l + Unh(p(p*))) (4.11) 

when 


- £<■ 


a cos (y 1 + V) j 
acot(y 1 + *y, 1 ) 


(4.12) 


for y* J < y* < oo (4.13) 


In equations 4.11 to 4.13, y c 2 is the half velocity point, 
£ in the length of the minor axis 8& it the momentum 
thickness on the minor axis and jn 2 is a real value 
satisfying: 


‘•nMffte*)) * 1 (414) 

For the subsonic low Reynolds number jet, Ho and 
Gumark, ref. 15, found very good agreement between 
the profile defined by equations 4.11 to 4.14 and ex- 
periment. Comparing figures 11a ar.d 13 shows that 
the velocity profile describes qualitatively the velocity 
curves in the physical jet cross section. 

For the present supersonic elliptic jet, the ofT axis 
mean velocity measurements, figure 11a, suggest that 
the mean velocity may be approximated accurately 
using a function which depends on one coordinate di- 
rection. While the experimental data was not used 
to generate the mean velocity profiles needed by the 
computations, it was used to approximate the cross 
sectional geometry of the jet. It is believed that data 
without the minor axis distortion should be used to 
approximate smoothly the mean velocity profile for a 
generalised supersonic shock-free elliptic jet. 

4.2 Numerical Analysis 

A hybrid numerical scheme is used to solve the sta- 
bility problem defined by equations 4.5 to 4.7. Hybrid 
techniques are numerical methods which combine se- 
ries approximations with finite difference calculations. 
The present scheme uses a pseudospectal series based 
on the Chebyshev polynomials to generate the dis- 
crete representation of the Rayleigh problem. 

In the present study, the pseudospectra] series ap- 
proximates the large scale pressure fluctuation, p, in 
the asimuthal direction, while the finite difference 
technique approximates p in the radial direction. On 
the edges of the jet cross section, the Rayleigh equa- 
tion reduces to the Helmholts equation. This equa- 
tion has exact analytic solutions in terms of Hankel 
functions of the first and second kinds. These solu- 
tions are used to define the boundary conditions on 
the edges of the jet. 

Moreover, the boundary conditions cutting through 
the jet, are based on the periodicity of p. It may be 


shown that either p or dp/ 88 is sero on the major and 
minor axes of the jet cross section. This allows the 
computational region to be restricted to one quadrant 
of the shea 7 jrer, as shown in figure 13. 

Once the boundary conditions cutting through the 
jet are applied to the pseudospectral series, the exact 
solutions on the edges of the shear layer are used as 
initial conditions to integrate the Rayleigh equation . 
The integrated solutions of the stability problem are 
then matched at the center of the shear layer using 
a generalised shooting technique. The matching pro- 
cess generates a matrix of solution vectors which are 
implicit functions of a fixed real frequency and some 
guessed complex wavenumber. This implies that the 
unknown wavenumbers, a e , may be determined by 
requiring that the determinant of the solution matrix 
be sero. The present method applies the iterative 
Newton-Raphson scheme to compute the wavenum- 
bers. 

Details of the hybrid numerical scheme and exam- 
ples of the stability computations for incompressible 
non-circular jets are given in ref. 3. 

5.0 COMPUTATIONAL RESULTS 

In this section the Rayleigh problem governing the 
linear inviscid isothermal stability of the compress- 
ible elliptic jet is solved using the hybrid numerical 
scheme described in the section 4.3. Calculations 
are performed to compute the wavenumbers for the 
jet at the axial stations located at X — 0.25L C and 
X = 0.50L C . These calculations sure performed for 
the four asimuthal normal modes corresponding to 
the varicose, flapping, and the ofT axis instabilities. 
The instantaneous pressure for these instabilities is 
shown in figure 14, ref. 14. In addition, the Strouhal 
numbers associated with the neutral stability point 
for the varicose instability are computed for the ax- 
ial stations between X = 0.25L C and X = 0.60I C at 
increments c* 0.05L C . 

5.1 Stability Analysis 

For all the numerical experiments, the aspect ratio 
of the half velocity point and the equivalent radius, 
defined by equation 4.5, have been fixed at: 

\ = 1-88 (5.1) 

and 

L = 1.0 (5.2) 

respectively. The value of 1.88 for the aspect ratio 
of the jet was choosen by plotting the aspect ratio of 
the half velocity point as a function of axial distance, 
figure 15. This aspect ratio describes an average shear 
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layer geometry for the jet which includes the effect 
of the bulge on the minor axis. It should be noted, 
however, that if the nossle imperfection were removed 
the aspect ratio for the jet would increase to a value 
■lightly greater than 2.0. 

In addition to fixing the aspect ratio and the refer- 
ence length scale, the momentum thickness must be 
determined for the calculations. From figure 12, it is 
seen that the momentum thickness is nearly constant 
on the major and minor axes down to the end of the 
potential core. Theretore, the momentum thickness 
for the minor axis is taken to be a linear function of 
axial distance. Then from figure 12, the momentum 
thickness for the calculations is defined to be 

6 b = 0.112(X/I e ) + 0.004 (5.3) 

for 0.25 < ( X/L e ) < 0.60. 

The program to solve the Rayleigh problem was ver- 
ified for incompressible and compressible examples in 
Baty, ref. 12. For all the jet cross sections considered, 
the hydrid method was run with 7 interior grid points 
and a finite difference step site on the order of 0.005. 
Figure 16 shows the variation of the axial growth rate 
as a function of frequency for the varicose instabil- 
ity of the jet at X = 0.25 L c . The variation of the 
phase velocity, given by u/Re(a t ), is shown in figure 
17. These results are typical of all the calculations for 
the four instabilities of the elliptic jet. To determine 
the most unstable modes for the different instabilities 
the three largest growth rates were interpolated using 
a second order polynomial. Table 1 shows the fre- 
quencies of the maximum growth rates for the elliptic 
jet for several instabilities at the two axial stations of 
X = 0.25l t and X = 0.50I e . 


lc 

Mode 

Freq 

Wavenumber 

0.25 

V 

4.00 

7.40 - 2.44s 

0.50 

V 

2.24 

3.85 - 1.71s 

0.25 

Fa 

3.90 

7.30 - 2.65s 

0.50 

Fa 

2.51 

4.24 • 1.96s 

0.25 

Fb 

3.86 

7.02 - 2.41s 

0.25 

Of 

4.41 

7.45 - 1.75s 


Table 1: Examples of the maximum growth rates and 
their associated frequencies for the Mach 1.52 elliptic 
are shown for the axial stations of X = 0.25L e and 
X — 0.501c. Here, V is the varicose instability; Fa 
is the flapping instability about the major axis; Fb is 
the flapping, instability about the minor axis; and Of 
is the off axis, peanut shaped instability. 

In addition to computing the most unstable modes, 
the neutral frequency as a function of axial length was 
predicted for the varicose instability. The varicose in- 
stability was choosen for the calcuation, because the 
experimental measurements in section 3.2 suggested 


that this was the main instability for the shock free 
elliptic jet. Figure 16 shows the neutral frequency for 
this instability at the axial station of X = 0.25£ c . 
As discussed in Liou, ref. 16, the neutral frequency 
at a given jet cross section represents the peak fre- 
quency of the large scale structure at that location in 
the jet. The only peak frequency measured experi- 
mentally was the frequency associated with the noise 
directed to the narrow cone along the jet axis from 
a region near the end of the potential core. This fre- 
quency corresponds to an approximate Strouhal num- 
ber of 0.19. Results from studies on axisymmetric jets, 
ref. 9, show that this frequency is close to the large 
scale instability frequency near the end of the poten- 
tial core. In figure 16, the predicted neutral frequen- 
cies are plotted as Strouhal numbers for axial length 
nondimensionalised with respect to the potential core 
length. 

5.2 Discussion 

Both the experiment and linear inviscid spatial sta- 
bility theory show that only the lowest order az- 
imuthal modes govern the stability of the supersonic 
elliptic jet. Experimentally, the varicose mode is 
found to dominate all other modes, particularly those 
frequencies that dominate sound radiation from the 
jet. Seiner, et al., ref. 17, shewed that for the round 
jet the varicose (i.e. axisymmetric) instability domi- 
nates all other instabilities for shock free jet plumes. 
Plumes that contain shocks are often dominated by 
flapping about a major axis. The elliptic jet results 
appear consistent with the previous round jet data. 

The experimental results in reference 9 and those of 
Troutt and McLaughlin, ref. 18, indicate that the fre- 
quency associated with peak amplitude emission for 
round shock free supersonic jets is assoicated with 
a region near the end of the potential core. In the 
present experiment, this measured frequency corre- 
sponds to a Strouhal number of 0.19. The analysis in- 
dicates, that for the varicose instability, this Strouhal 
number dominates the flow in a region near 60 per- 
cent of the potential core. This feature of the linear 
stability analysis has been observed previously, ref. 
9, for round jets, where it is conjectured that large 
scale waves grow more rapidly in an initially laminar 
jet column. The spatial stability analysis, which as- 
sumes a laminar shear layer with a superimposed large 
scale fixed frequency wave, has always achieved better 
agreement to experimental results for lower Reynolds 
number jets. What effect the simulation of mean ve- 
locity profiles has on the spatial stability analysis is 
unclear, except as noted in reference 9, the analysis is 
extremely sensitive to second order derivative proper- 
ties of the assumed mean velocity profile. 

The substantial growth in the major axis momen- 
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turn thickness at the end of the potential core is pos- 
sibly related to decay of the most amplified wave in 
the jet column. An unfortunate consequence of the 
present is that it cannot predict such behavior. There 
is a strong need to understand the mechanism behind 
such behavior, since control of the jet’s spread rate 
could possibly be affected through manipulation of 
the small nonlinearity controlling the major axis mo- 
mentum thickness. 

0.0 CONCLUSIONS 

This paper presents the results of linear inviscid 
spatial stability theory used to calculate the most 
highly amplified modes of the large scale coherent 
structures in the supersonic shock free elliptic jet. The 
results of the analysis have shown that only the lowest 
order instabilities are significant and that the frequen- 
cies for the maximum growth rate decrease with axial 
distance from the nossle. In a parallel effort an ellip- 
tic nost le with an exit aspect ratio of 2 was studied to 
verify the analytical results. This nossle was designed 
for an exit Mach number of 1.5. The experimental re- 
sult* indicated that the varicose instability was the 
dominant «'<atiai mode when operating the nossle at 
its shock free point, a result consistent with round 
jet observations. This result is in agreement with the 
theoretical analysis in this paper. The measured fre- 
quency associated with peak amplitude noise radia- 
tion is found to agree reasonably with that frequency 
predicted to dominate the large scale structure near 
the end of the potential core. 

The present results provide encouragement for use 
of the inviscid Rayleigh equation to provide a rea- 
sonable model to account for dynamical processes 
associated with non-circular jets. To improve the 
present quantitative agreement with experiment may 
require a better description of the real mean flow. In 
its present form, the determination of a fixed fre- 
quency wave’s neutral point is computationally in- 
tensive. This suggests that more efficient numerical 
methods are needed to solve the spatial stability prob- 
lem. 

The authors wish to thank Mr. Bernie Jansen for 
his help in acquiring the aerodynamic data associated 
with this research. The authors also wish to thank Dr. 
Glenn Laguna of Sandia National Labs for his help in 
setting up the batch processes needed to perform the 
compntational experiments. 
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8.0 FIGURES 
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Figure 1 Internal geometry of the elliptic nostle 
major and minor axes cross sections 



Figure 2a 



( igure 2b 


Figures 2a-b Schlieren photos shewing the major 
and minor axes of the elliptic jet operating at the 
design conditions 
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Figure S: Elliptic near field microphone array. 



40 SS » 25 20 IS 10 


R/Daq 


• danota* microphona 




Figure 5: Contour map of OASPL data, M t — l.S, 
Mj = 1.52. 
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Figure 4: Seven radial portion* for the linear micro- 
phone array. 


Figure 6: Example narrow band acoustic spectrum 
in major noise producing region along minor axis. 
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Figure tc 

Figure 7: Noiile exit seer field pressure apeetram in 
minor axis plane. 
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Figure h 
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Figures 8a*c: Asimnthal distribution of 1/S octave 
band crow correlation coefficients with modal ampli- 
tudes: a. 1000 Hs band, b. 16S0 Hs band, and c. 
2000 Hs band. 
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Figure 9: Centerline distribution of mean axial ve- 
locity. 


Figure 8b 
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Figure IS: Jet k id er o w nectioa at 0.25L* ahow- 
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Figures 14»4: The axunathal iastaateaeoas ptes- 
nut conditions for the famt iastabilities: a. vancoae, 
b. flapping aboat the major axis, e. off axis (peaaat 
shaped), aad d. flapping aboat the aiiaor axis. 


figare IS: Aspect ratio of the elliptic jet's half ve- 
locity point as a fuctioa of axial distaace. 
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Figure 16: Variation of the axial growth rate with 
frequency. Varicose instability for the elliptic jet at 
X ~ 0.25U C 
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Figure It: Variation of the centra] Stroahal Bomber 
with axial length for the varicose instability 




